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Abstract. We have carried out a study of the orthogonal polarisation mode behaviour as a function of frequency of 18 pulsars, 
using average pulsar data from the European Pulsar Network (EPN). Assuming that the radiation consists of two 100% polarised 
completely orthogonal superposed modes we separated these modes, resulting in average pulse profiles of each mode at multiple 
frequencies for each pulsar. Furthermore, we studied the frequency dependence of the relative intensity of these modes. We 
found in many pulsars that the average pulse profiles of the two modes differ in their dependence on frequency. In particular, 
we found that pulse components that are dominated by one mode tend to increase in intensity with increasing frequency with 
respect to the rest of the profile. 
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1. Introduction 

Single pulse studies of the position angle (PA) of linearly po¬ 
larised radiation from a number of pulsars show that it is built 
up of two modes of polarisation, which are s e parated in a ngle 

jitinebring et alJ 

1199211 . It is believed 
that these orthogonally polarised modes (OPM) reflect the 
eigenmodes of the magneto-active plasma in the open mag¬ 
netic field-lines above the pulsar polar cap. Three modes of 
wave propagation are allowed in this region, one of which is 
the ordinary sub-luminous mode, which cannot escape the pul¬ 
sar magnetosphere due to Landau damping and is thus of no 
interest to the present work. The remaining two modes are the 
ordinary super-luminous mode (O-mode) which is polarised in 
the local plane defined by the external magnetic field and the 
wave-vector, and the extraordinary mode (X-mode), which is 
polarised perpendicular to this plane IXrons & Barn aril 98(4 
IPetrovalkoO f;). According to iBamard&AronTlT^Hfil) refrac- 
tion can separate the X- and O-modes beams, which have dif¬ 
ferent indexes of refraction, by many beam-widths. The ten¬ 
dency for average profiles to have constant widths above a crit¬ 
ical frequency is then explained by this separation occurring 
above a critical height dSieber et alii 1975I) . Two conditions are 
required for this: (1) the radio emission mechanism has to be 
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broadband in frequency over a narrow range of heights above 
the stellar surface and (2) the gradient in the plas ma den sity 
has a transverse component to the radial direction IMcKinnonl 
I1997I) . The independent propagation of the two modes in the 
open flux zone can produce the abrupt orthogonal transitions in 
polarisation position angle that are commonly observed in stud¬ 
ies of individual pulse polarisation ( Manchester et alJll975l) . 
The merging of the beams might also account for depolar¬ 
isation o f pulse average profiles with increasing radio fre - 
quency dMcKinnonl 1997tIvon Hoensbroech~& Xilourisll997l) . 
Alternatively, the OPM transitions might not be due to refrac¬ 
tion, but due to switching between significant and insignificant 
conversion of O-mode into X-mode (lPetroval200ll) . 

One of the most important questions about OPM is whether 
the polarisation modes are disjoint or superposed. For disjoint 
modes the polarisation of the observed radiation is given by ei¬ 
ther one of the two modes at each point in time. Superposed 
modes occur simultaneously and the polarisation of the ob¬ 
served electromagnetic radiation is given by the vectorial ad¬ 
dition of the Stokes parameters Q, U and V of both modes. 
fCordes et all d 1978l) . who first posed this question based on the 
polarisation from PSR B2020+28, assumed the modes to be 
disjoint and stated three arguments for this assumption: (1) the 
degree of polarisation is fairly steady from one pulse to the 
next, (2) if the modes are superposed, one might expect oc¬ 
casional sign changes in the complex value Q + iU when the 
linear polarisation is low, due to random fluctuations in the rel- 
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ative strength of the two modes and (3) the instantaneous mode 
of polarisation seems correlated with the total intensity, which 
inconsistent with disjoint occurrence of the modes. However, 
IStinebring et al.l f 1984) found that the modes appear to be su¬ 
perposed by studying sensitive polarisation observations of a 
number of pulsars. It was found that average values of the 
fractional linear polarisation were small at longitudes where 
both polarisation modes occurred with nearly equal frequency. 
The same effect was later found in the instantaneous values of 
the fractional linear polarisation in single pulse observations 
of PSR B2020+28 JMcKinnon & Stinebrindlil998h . This phe¬ 
nomenon is present on timescales of hundreds of microseconds 
to hours, which is what one would expect for superposed OPM. 
For disjoint modes one would expect the fractional polarisation 
to be either high or low, depending on the fractional polarisa¬ 
tion of the active mode. 

Assuming that the polarisation modes are superposed, that 
they are completely polarised and completely orthogonal, the 
average intensity of the individual polarisation modes as a func¬ 
tion of longitude can be determined from the Stokes parameters 
of the average pulsar signal (IMcKinnon & Stinebrina200d) . In 
this paper we make these assumptions and determine the indi¬ 
vidual polarisation profiles of 18 pulsars as a function of both 
longitude and frequency (using data from the EPN database). 
In section 13 we give a description of our analysis, the results 
of which are shown in section |3 In section 0 we discuss the 
results and finally we give our conclusions in section|3 

2. Method 

McKinnon and Stinebring have developed a method to sep¬ 
arate the modes of electromagnetic radiation that is built up 
out of two 100% polarised completely orthogonal modes that 
are superposed IMcKinnon & Stinebringl ; 200Q ). This method 
is illustrated by means of vectors in Fig. (TJ Each pulse longi¬ 
tude will have a value for Stokes Q, U and V, which can be 
expressed by a polarisation vector in Q-U-V-space (Poincare 
space). The intensity corresponds to the radius of the sphere. 
Polarisation vectors lying on the sphere indicate full polarisa¬ 
tion, whereas vectors lying within the sphere indicate partial 
polarisation. The fraction of polarisation is given by the length 
of the polarisation vector divided by the radius of the sphere 
(the intensity). Furthermore, two vectors with orthogonal po¬ 
larisation will be antiparallel. The method first fixes the orien¬ 
tation of the polarisation vectors to lie in the Q-V plane, shown 
in Fig.Q by rotating the vector around the V axis. In this figure, 
the average intensity of the observed radiation is given by the 
circle with radius I e , the observed polarisation vector is given 
by P e and the underlying polarisation of the orthogonal modes 
are given by P\ and Pi. There are now two relationships be¬ 
tween the observed intensity and polarisation vector ( I e and 
P e ) and the polarisation vectors of the orthogonal modes (Pi 
and Pi). First, I,, is equal to the sum of the lengths of Pi and Pi 
and second, P e is equal to the vectorial sum of Pi and P 2 . With 
these relationships the two orthogonal polarisation vectors can 
easily be determined from the Stokes parameters. 

We apply this method to 18 pulsars for which polarimetric 
average pulse profiles are available in the EPN database. Given 



Fig. 1 . Poincare circle that illustrates the method used to de¬ 
termine the intensity and polarisation of two orthogonally po¬ 
larised modes from the observed radiation. Pi and Pi repre¬ 
sent the antiparallel polarisation vectors of both modes. The 
radius of the circle, I e , represents the total intensity of the ra¬ 
diation, which is equal to the sum of the lengths of vectors P\ 
and P 2 , because the polarisation modes are superposed inco¬ 
herently. P e represents the polarisation of the electromagnetic 
radiation which is equal to the vectorial sum of P\ and P 2 . Note 
that the polarisation vector of the observed radiation has been 
rotated around the V axis to lie in the Q-V plane. 

the assumptions that the radiation is built out of two superposed 
orthogonal modes that are 100% polarised, the separation is 
completely correct. Because the superposition is linear, sepa¬ 
ration of the average polarisation profile at each longitude will 
result in the means of the individual modes at that longitude. 
Because of this, we do not see the polarimetric behavior of sin¬ 
gle pulses, nor any polarimetric behavior on small timescales. 
If, for example, the separation from the average polarisation 
profile at a fixed longitude would reveal both modes to be of 
equal strength, then one solution might be that at this longi¬ 
tude all single pulses are built out of the superposition of two 
modes with equal strength. Yet it might also be the case that 
at this longitude one half of the single pulses would be built 
out of just one of the orthogonal modes and the rest of the sin¬ 
gle pulses would be built out of just the other mode. The data 
from the EPN database consists of time-averaged data at multi¬ 
ple frequencies. Of these pulsars, 11 are known from previous 
studies to show OPM. It should be noted that this method gives 
separated intensities corresponding to both modes at each lon¬ 
gitude, but there is no direct knowledge as to determine which 
value belongs to which mode. We therefore have to be cautious 
when making a longitude plot of the intensity of both modes 
as there are always two possibilities when connecting intensity 
points from one longitude to the next. Generally, it is safe to 
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assign values with the highest intensity to one mode and val¬ 
ues with the lowest intensity to the other. When these values 
lie very close to each other, the assignment of value to modes 
is made on the basis of visual inspection for changes in modal 
dominance as evidenced by sharp jumps of nearly 90° in the 
PA. In some cases the intensity points are connected by finding 
a smooth line through the values. We have defined the polari¬ 
sation mode with the highest integrated intensity at the lowest 
available frequency as the strongest mode. We thus produce the 
intensity of the separated modes, as well as the PA and elliptic- 
ity angle (^atan(V/ ^Q 2 + U 2 ) as a function of frequency. 

Furthermore, we have calculated the ratio of the longitude- 
averaged intensity of both modes at each frequency. For pul¬ 
sars with two pulse components this has been done for both 
components separately and for both components added to¬ 
gether. The error on the average intensity was calculated as 
crj = ~yN ■ cr, where N is the number of bins of the pro¬ 
file and cr is the noise level, given by the standard deviation 
of the noise as estimated from data in the off-pulse region of 
the profile. For 14 pulsars, the pulse profiles at all frequen¬ 
cies were aligned by visual comparison of the profiles with 
the corresponding plots from iKuzmin et al.1 (1 998 ). When a 
profile at a particular frequency had no direct comparison, 
the alignment was estimated from the frequency development 
from the plots from lKuzmin et al.l ( 1998 1. Pulsars B0144+59, 
B1039-19, B1800-21 and B1944+17 were not listed in this 
paper and were aligned by comparing pulse components and 
characteristics in the polarisation from frequency to frequency. 
We can therefore not guarantee that these pulsars are aligned 
correctly. 

3. Results 


L 



Line of Sight 


Line of Sight 


Fig. 2. Schematics of two field lines from which radio waves of 
two different frequencies are observed. L is the line that con¬ 
nects the locations where the dipolar field lines are directed to¬ 
wards the observer. Due to radius-to-frequency mapping, low 
frequency radio waves are emitted higher in the magnetosphere 
than high frequency radio waves. As the image shows, the low 
frequency radio waves that can be observed are emitted on 
field lines that intersect the surface closer to the magnetic axis 
than the field lines that emit the observed high frequency radio 
waves. 


If we assume that emission altitude is inversely related to the 
frequency of emission (radius to frequency mapping) and that 
the magnetic field-lines are dipolar shaped at the emission 
height, observing at different frequencies has two implications: 
(1) radiation at high frequency arises lower in the pulsar mag¬ 
netosphere than radiation at low frequency and (2) radiation ob¬ 
served at high frequency is emitted at field-lines that intersect 
the surface further away from the magnetic axis than radiation 
observed at the same pulsar phase at low frequency (see Fig.E}- 
When studying the average profiles at multiple frequencies, this 
should be kept in mind. 

Figures [4] through eh show the results of the mode sep¬ 
aration at all frequencies for all pulsars. The left columns 
show the average total intensity (solid line) and the average 
intensity of both modes (dotted and dashed lines), the mid¬ 
dle columns show the PA of the average Stokes parameters 
and the right columns show the ellipticity angle at all frequen¬ 
cies. The PA and ellipticity angle at pulse longitudes where 
the polarised intensity is less than twice the RMS noise are 
not shown. Figures E3 through |M] show the frequency devel¬ 
opment of the ratio of the average intensities of both modes, 
plotted on a log-log scale, for all pulsars. The straight line 
through the points was determined by the least squares fitting 
technique. Whenever these points show little deviation from 
the straight line fit, the related power law index is shown in 


Table El As there is no preferred mode, the values are always 
taken to be positive. The point-style indicates the reference cor¬ 
responding to the source of the data in the EPN, which is listed 
below Fig. El The polarisa t ion of the 1.41 GHz data from 
Ivon Hoensbroech & Xilouris ill 9971) was frequently found to 
differ from other 1.41 GHz EPN data and from the develop¬ 
ment of the polarisation as a function of frequency. Thus these 
data were not included in the analysis. 

Let us now consider each pulsar in detail. For some pul sars 
we will refer to the classification by Ivon Hoensbroech et all 
dl998h . who have classified a number of pulsars into four 
groups, named after one of the pulsars in that group. The first 
group of pulsars is exemplified by type 0355+54, for which 
one component is highly polarised, whereas the rest of the pro¬ 
file has little polarisation. The component that is highly po¬ 
larised has a flatter spectral behaviour than the rest of the pro¬ 
file and therefore dominates at higher frequency. The second 
group consists of 0525+21-like pulsars, which, in the classifi¬ 
cation system from Rankin ( 1983 ). are all examples of conal 
doubles. They have two moderately linearly polarised com¬ 
ponent the PA follows the rotating vector model (RVM) of 
iRadhakrishnan & Cooke ( i 19691 nearly perfectly and there is 
usually a core-type bridge between the components. The third 
group consists of 1800-21-like pulsars, which are young, have 
a very high spin down power, E, a relatively flat radio spectrum 
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behaviour, detectable X-ray emission and radio emission over a 
wide fraction of the pulse period. They’re also nearly fully po¬ 
larised, linearly and circularly and do not show the commonly 
observed effect of depolarisation towards high frequencies. The 
final group consists of 0144+59-like pulsars, which show an 
unusual polarisation behaviour. Their degree of circular polari¬ 
sation increases strongly with frequency. 

B0144+59 

Pulsar B0144+59 defines the class of 0144+59-like pulsars, 
as described above. Figure 0 shows that at all frequencies ex¬ 
cept 4.85 GHz there is a mode-switch around -2°, which is 
associated with a jump in the PA. At 0.4 GHz the first com¬ 
ponent can be associated with one mode, whereas the second 
and third components seem to be associated with the other 
mode. Although the intensity fluctuations around 2°, which 
suggest a second and third component, are not significant and 
might be just one component. The association of the modes 
with components is also valid at higher frequencies, even when 
the first component becomes weaker and almost disappears at 
4.85 GHz. At this frequency the remaining component is al¬ 
most completely dominated by one mode. Also, the slope of the 
position angle traverse increases. This pulsar shows an exam¬ 
ple of a component (left component) being dominated by a sin¬ 
gle mode and slowly disappearing with increasing frequency. 
It would be interesting to see whether this trend continues to¬ 
wards even lower frequencies. Unfortunately, such data was not 
at our disposal. The frequency development of these modes is 
nicely fit by a power law. 

B0355+54 

Pulsar B0355+54 defines the class of 0355+54-like pulsars, as 
mentioned above. Figure 0 shows that the leading component 
consists of mainly one mode (the stronger mode) and the sec¬ 
ond component is dominated by the other mode. At all frequen¬ 
cies, except perhaps for 10.55 GHz, the PA jumps and the el- 
lipticity angle switches sign at the same longitude where the 
modes change dominance. At low frequency the dominance of 
the weaker mode is only a few degrees in longitude, but be¬ 
comes broader with increasing frequency. This mode can be 
associated with the second component, but the stronger mode 
is present there as well. The first component rises in inten¬ 
sity at higher frequencies, relative to the intensity of the first 
component, which contains only one mode. This behavior has 
already been noted bv lvon Hoensbroech et al. ( 1998 1 The PA 
sweep of the left component is straight and has a small slope, 
which slightly increases with increasing frequency. The right 
component shows the 90° jump, after which the sweep remains 
straight for a few degrees (the number of degrees depending on 
frequency) and then has a smaller jump, followed by a curved 
sweep. The second jump in the PA can be associated with a 
second change in dominance of the two modes. However, at 
high frequencies this is not clear. It is also not clear whether 
the curve at the trailing edge of the PA is the natural PA sweep 
of one of the modes, or that the curve is the result of a (non- 


orthogonal) superposition of both modes. The frequency devel¬ 
opment of the modes is nicely fit by a power law, except for a 
jump at 1.642 GHz. 

B0450+55 

Pulsar B0450+55 falls into the same group of pulsars as 
B0355+54 Ivon Hoensbroech et alll998l) . Figure[6]shows that 
at a frequency of 0.910GHz the PA jumps 90° twice around 
0° longitude, indicating a brief mode-change which is visible 
in the average profile of the modes. The first jump shows a 
gradual change. The ellipticity angle switches sign at the same 
longitude in accordance with the mode-change. At lower fre¬ 
quencies the PA also jumps twice at the same longitudes, but 
this does not show up in the average profile of the modes. These 
jumps are less than 90° and also show a gradual change. It is 
interesting that at 0.404 GHz the first jump is upwards while at 
0.610GHz the first jump is downwards, while the rest of the 
PA sweep is very alike. The first component of the average in¬ 
tensity is strongly dominated by one mode and the intensity of 
this component increases with frequency relative to the second 
component, in which both modes are present and nearly equally 
strong. Because it is dominated by one mode, Fig.[32]does not 
show the frequency development of the modes in the left com¬ 
ponent. For both the right component and the entire profile, the 
frequency development of the modes do not seem to be fit by a 
power law. 

B0525+21 

Pulsar B0525+21 defines the class of B0525+21-like pulsars, 
as described above. Figure 0s hows that both modes are active 
in both components. At high frequency, there is strong depolar¬ 
isation in both components, which is represented by the weaker 
mode becoming stronger at higher frequency. The modes in 
both components behave very similarly, although the ratio be¬ 
tween integrated intensity of the stronger and weaker mode is 
higher in the right component. The intensity of this compo¬ 
nent increases with frequency, relative to the left component. 
The PA sweep at 4.85 GHz shows a 90° phase jump, suggest¬ 
ing a mode-switch somewhere in the bridge between the two 
components. However, the frequency development of the PA 
sweep clarifies that this jump is actually the result of the nor¬ 
mal RVM. The ellipticity angle gives little information about 
mode-switching, because of the low circular polarisation. This 
leaves us with no indication of mode-switching in this pulsar. 
The frequency development of the modes are nicely fit by a 
power law in both components. 

B0809+74 

The polarisation of pulsar B0809+74 has been well studied in 
the past. Recently, [Rankin et~al . ( 2004 has shown that the two 
orthogonal modes of polarisation are the primary reason for 
a variation in the longitudinal separation between the drifting 
subpulses of pulsar B0809+74. Furthermore. lEdwardsi ( 2004 1 
has confirmed that the polarisation fluctuations are the result 
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of two orthogonally polarised modes, although an apparently 
randomly polarised component needs to be superposed. The 
pulsar has two components. The first component arises with 
increasing frequency and becomes visible around 0.925 GHz 
and higher and is almost completely dominated by the stronger 
mode (see Fig. E} In the second component both modes are 
almost equally strong. At frequencies of 0.606 and 1.408 GHz, 
the PA jumps by approximately 60° and 80°, respectively. This 
does not seem to be associated with a mode-change. There is 
no clear frequency development of the modes. 


B0823+26 

Pulsar B0823+26 has only one component, in which both 
modes are active (see Fig. [9}. Both the profile and the relative 
intensity of the modes do not change very much over frequency. 
Around the middle of the profile there is a small jump in the PA 
sweep, which increases in size with increasing frequency. At 
the same longitude as this jump the modes are almost of equal 
strength. Also, the weaker mode becomes stronger with respect 
to the stronger mode at higher frequency. The frequency devel¬ 
opment of the modes are nicely fit by a power law. 


B0834+06 

Pulsar B0834+06 consists of two components with a high 
bridge connecting them. From studies of individual-pulse po¬ 
larisation observations it has been concluded that both modes 
are equally strong in the two componcntsGiut only one mode 
dominates across the bridge Istinebring et all'19841) . However, 
Fig. EH shows that the difference in average intensity between 
the modes is small across the bridge as well. The PA sweep 
at 0.408 GHz and 1.408 GHz are very similar, however, the 
PA sweep at 0.610GHz has a drop in the center of the profile. 
Except for one point with a large error, the frequency develop¬ 
ment of the modes can be fit by a power law. But since there 
are only four points in total, this is not very significant. 


B0919+06 

The average profile of pulsar B0919+06 changes greatly over 
frequency. At low frequency there seem to be two components 
visible, as well as a slope at the leading edge. With increasing 
frequency the two components merge into one and the slope 
becomes smaller and even disappears at 4.85 GHz. The pro¬ 
file is dominated by one mode, with the exception of the slope 
and part of the left component, in which both modes are al¬ 
most equally strong (see Fig. ED- The weaker mode moves 
to the right with respect to the stronger mode, at rising fre¬ 
quency. This results in the extra component at low frequency 
(0.408 GHz), associated with the weaker mode, to disappear 
towards higher frequencies. There is no indication of mode¬ 
switching. The frequency development of the modes seem to 
be fit well by a power law, with the exception of the datapoint 
at 4.85 GHz. 


B0950+08 

Pulsar B0950+08 is remarkably similar to pulsar B0919+06. 
They both have one dominant component a weak second com¬ 
ponent at low frequencies and a slope at the leading edge. 
However, the polarisation modes of pulsar B0950+08 are al¬ 
most equally strong over the entire frequency range and do not 
show any drift (see Fig. H3>. There is no evidence of a mode- 
switch in the average profile. The PA does rise and fall where 
the modes have nearly equal power, but only by 70°. The rise in 
the PA becomes broader with increasing frequency and at 1.642 
and 4.850 GHz, the PA has changed significantly. At 1.642 GHz 
it is clear that the PA drops instead of rises and by a smaller 
amount than at lower frequency. The ellipticity angle falls and 
rises at the same longitude at which the PA changes. This dip 
also becomes broader towards higher frequency. There does not 
seem to be a clear frequency development of the modes. 


B1039-19 

Pulsar B1039-19 is B0525+21-like. Figure[j3]shows that both 
modes are active in both components and that the PA does in¬ 
deed follow the RVM. At higher frequency, the weaker mode 
becomes stronger. At low frequency the PA sweep shows a typ¬ 
ical RVM curve, but towards higher frequency the sweep starts 
with a rise and then drops in a straight line. For the first com¬ 
ponents, the frequency development of the modes seems to be 
fit by a power law. For the second component the fit is not so 
good, but the same trend is visible. 


B1133+16 


This pulsar has two components. In the left component both 
modes are almost equally strong (see Fig. E3- The right com¬ 
ponent is dominated by one mode at low frequencies, but the 
modes become almost equally strong at higher frequencies. 
This behavior is the opposite of what would be expected by 
comparing the results on the polarisation of B1133+16 from 
Backer_&_Ranki3] 198(1) (hereafter BR) with the results from 
Stiiiebring et alJ (119841) (hereafter SCRWB). It can be seen 
there, that at low frequency the right component shows OPM 
(Fig. 9 in BR), which becomes less visible at higher frequency 
(Fig. 17 in SCRWB). However, OPM are expected to occur 
when both modes are about equal in strength which, according 
to our results, is at high frequency rather than at low frequency. 
This could indicate that the method for separating the modes 
is not succesful in this case. However, we must keep in mind 
that we are only showing average profiles. If the (refracted) 
path of one of the polarisation modes in the magnetoactive 
plasma can vary in time, then the strength of the polarisation 
modes can briefly change, possibly making OPM visible in sin¬ 
gle pulse studies. At 1.408 GHz the PA jumps upward by 90° 
and jumps downward again one degree later, which results in 
a mode-change, as can be seen in the plot of the average in¬ 
tensity of the modes. At 1.642 GHz this jump is still present, 
but has become more gradual, is less than 90° and does not re¬ 
sult in a mode-change. At 1.71 and 4.85 GHz these jumps are 
also present, but are first downward, then upward. At 4.85 GHz 
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the average intensity profile of the modes show a brief mode- 
change. The modes in the left component seems to have a clear 
frequency development, whereas the modes in the right com¬ 
ponent do not. 

B1237+25 

Pulsar B1237+25 is known for its profile “mode-switching” 
(which means that it has two different average profiles, called 
“modes”, which are not to be confused with the orthogonal po¬ 
larisation modes) (Bart el et al.lll982() . By comparing th e EPN 
profiles of pulsar B1237+25 with those in (Bartel et al. 1982), 
we conclude that the abnormal mode can only be faintly present 
in the average pulse profiles. There are many components in 
this pulsar, which are all dominated by the stronger mode (see 
Fig. EJ. While the weak component at a longitude of -3° 
can be associated almost completely with the stronger mode, 
the weak component at a longitude of H5 can be associated 
with the weaker mode. These associations are strongest at low 
frequencies. At higher frequency, the weaker mode becomes 
stronger with respect to the stronger mode. According to BR, 
this pulsar shows OPM in the two main components at a fre¬ 
quency of 0.430 GHz, which, as in the right component of pul¬ 
sar B1133+16, is not consistent with the result of one domi¬ 
nant mode. At almost all frequencies, the position angle rises 
and falls along the bridge between the two dominant compo¬ 
nents. At this same longitude, there is a lot of circular polarisa- 
tion. This discontinuity in the PA is not associated with OPM 
IStinebring £Ldil5>&4 , but it does occur at a longitude where 
the weaker mode becomes stronger with respect to the stronger 
mode. Both the frequency development of the intensity of the 
modes in the left and right half of the profile follows a power 
law, but the modes in the right half have a very low power law 
index. 

B1737-30 

Pulsar B1737-30 has only one component and fallsjnto the 
sam e group of pulsars as B0144+59 Ivon Hoensbroech et al.1 
1 998 ). At a frequency of 1.408 GHz and higher, the component 
is almost completely dominated by one polarisation mode. At 
the trailing edge of the component the amount of circular polar¬ 
isation increases. This behaviour is intrinsic to the modes them¬ 
selves, as there is mainly one mode active throughout the pro¬ 
files. The frequency development of the intensity of the modes 
seem to follow a power law. 

B1800-21 

Pulsar B1800-21 defines the class of 1800-21-like pulsars, as 
described above. At low frequency this pulsar has up to 4 com¬ 
ponents, but at high frequency this pulsar has only two com¬ 
ponents, both dominated by the stronger mode (see Fig. E3 
The first component arises at a frequency of 1.408 GHz. At 
1.408, 1.56, 1.642 and 4.852 GHz there is a gap in the PA due 
to the low signal in the bridge between the two main compo¬ 
nents. The PA at 4.85 GHz shows the PA sweep over this bridge 


and clarifies that the stronger mode is dominant in both com¬ 
ponents. There is, however, a brief mode-switch in the bridge 
between the two main components at 1.408 GHz and another 
one at the trailing tail of the second component. At frequen¬ 
cies below 4.85 GHz the peak of the weaker mode is separated 
by approximately 15° (= 5.6 ms) from the peak of the stronger 
mode. At a frequency of 4.85 GHz, the weaker mode has almost 
completely disappeared. There is no clear frequency develop¬ 
ment of the intensity of the modes, but it is clear that the weaker 
mode becomes weaker with increasing frequency. 

B1929+10 

The profile of pulsar B1929+10 changes greatly over fre¬ 
quency, but at all frequencies there is only one dominant mode 
(see Fig. I18> . There is OP M at th e trailing part of the main 
component (Stinebring et al .111 9841) . which coincides with the 
rise of the weaker component. At 4.85 GHz the position angle 
suggests that there is a mode-switch in the leading edge of the 
profile. Due to the low intensity it is difficult to confirm this 
with the intensity plots. It is interesting that the weaker mode 
retains its shape and position, whereas the components of the 
stronger mode change greatly in shape and position. This is es¬ 
pecially evident at 4.85 GHz where one might naively expect 
the two components in the average intensity profile to lie at the 
same longitude as the second and third components of the 1.71- 
GHz profile, rather than to coincide with the first and second 
components. However, using the alignment from lKuzmin et all 
( 1998 ) the weaker mode does not change in longitude over fre¬ 
quency, which suggests the alignment to be correct. There is no 
clear frequency development of the intensity of the modes, but 
it seems that the weaker mode becomes stronger with increas¬ 
ing frequency. 

B1944+17 

Pulsar B1944+17 switches modes at the leading and trailing 
edge of the pulse at 0.610 and 1.408 GHz (see Fig. ll9t . Overall, 
the weaker mode has different intensity at different frequencies. 
It is roughly 1/4 of the dominant mode at frequencies of 0.410 
and 0.925 GHz and about 1/2 of the strength of the dominant 
mode at frequencies of 0.61, 1.408 and 4.85 GHz. However, 
there is no clear frequency development of the intensity of the 
modes. 

B2016+28 

The profile of pulsar B2016+28 has two components, which 
become less distinct towards higher frequency (see also 
dlzvekova et al.[ll993 l)). Fig. [3UI shows that these components 
can be associated with the two different modes. At low fre¬ 
quency, the weaker mode consists of two components and 
of only one component at higher frequency. The peak of the 
modes move towards each other with increasing frequency and 
become of almost equal strength. At frequencies higher than 
1.408 GHz, the shapes of the two modes are almost similar. Due 
to this behavior, the left part of the profile is dominated by one 
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mode at low frequencies while it is composed of two equally 
strong modes at hi gh frequency. This pu l sar is known for 
its m ode-changing iStinebring et al . 1984: [Arons & Barnardl 
Il986h . This can be seen at a frequency of 0.408 GHz, where the 
modes switch dominance. At the same longitude the PA jumps 
by approximately 90°, however the PA at 1.408 GHz has a sim¬ 
ilar behavior and does not seem to be associated with a mode- 
change. This suggests that the jump in the PA at 0.408 GHz and 
1.408 GHz are due to the RVM and respresent the PA sweep of 
the dominant mode. Moreover, at 0.610 and 0.925 GHz the PA 
jumps downwards and does not seem to be associated with a 
mode-change. It therefore remains unclear what the true PA 
sweep of this pulsar is. There does not seem to be a frequency 
development of the intensity of the modes. 

B2020+28 

Pulsar B202 0+28 is known for its abrupt mode-change 
i Cordes et alJ ll 9781 IStinebring et al.l 179841 lArons & Barnard 
Il986l) . This can be seen in the PA at all frequencies (see 
Fiu. 12 1k This pulsar also shows a large drift of the modes with 
respect to each other. Fig 0 shows the difference in longitude 
between the peaks of the two modes as a function of frequency. 
(The large errors are because of the uncertainty in the position 
of the peak of each mode.) The line through the points is a fit of 
the function a log(/) +b of the points, where / is the frequency 
in GHz and a and b are fit parameters. The values of these fit 
parameters are a = CF56±(T11 and b = 8?54±0? 13. There does 
not seem to be a frequency development of the intensity of the 
modes. 


4. Discussion 


Data from the EPN is very suitable to apply the mode- 
separation technique for superp osed modes, as prescribed by 
McKinnon & Stinebring ( 2000 ). This makes a frequency anal¬ 
ysis of OPM possible. Applying this technique to EPN data has 
three drawbacks (1) the pulse-profiles at different frequencies 
are not aligned, (2) the information on the observed pulsar flux 
is seldom available and (3) the data only show average pro¬ 
files, which might not represent the behavior of single pulses. 
Still, the results described here, give a clear picture of how the 
average properties of the two modes behave as a function of 
frequency. 

In Table Q we summarise the results of the analysis of the 
orthogonally polarised modes of the 18 pulsars, as well as some 
of the known parameters. W 50 and w 10 are the pulse widths at 
50% and 10% of the peak flux, respectively, S 400 and 5 1400 are 
the flux densities at 0.4 and 1.4 GHz respectively and E is the 
spin-down luminosity, given by -700, where 7 is the pulsar’s 
moment of inertia for which we take 10 45 g cm 2 and Q and Cl 
are the rotation frequency and the time derivative of the rota¬ 
tion frequency of the pulsar, respectively. Values for the period , 
pulse widths and flux densities are fro m [TavTo r et al.l ( 1993 ). 
Type and a and /j parameters are from iRankint l i 19931) . except 
for pulsars B0144+59 and B1800-21 for which we have deter¬ 
mined the type ourselves from the polarisation characteristics. 
The types are based on the classification system from lRankinl 
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Fig. 3. Separation in longitude between the peaks of the aver¬ 
age intensity of the two orthogonally polarised modes of PSR 
B2020+28 as a function of frequency. The dashed line is the 
best straight line fit through the points. 

il983t) which assumes a core radiation zone, surrounded by a 
conal radiation zone. They have the following meaning. S,/ is 
a single component profile; the line of sight grazes the outer 
cone of the polar cap. D is a double component profile; the line 
of sight cuts through the outer cone twice, but misses the core 
component. T is a triple component profile; the line of sight 
cuts through the outer cone twice as well as the core. S, is a 
single component profile; in this case the outer cone is miss¬ 
ing and the line of sight only cuts through the core. Finally 
there are two types that, in this view, can only be explained by 
two conal radiation zones around the core component; they are 
M, which is a multiple component profile (containing usually 5 
components) and T c , which is a triple component profile, which 
consist of three conal components. 

Let us now consider some of the observed phenomena. 

Mode-changing 

In almost half of the pulsars, there is a jump of 90° in the PA, 
which can always be associated with a change in the dominance 
of the modes. (See first row of Table Q. This can be during a 
fraction of the pulse period, as in pulsar B1133+16, or for most 
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Table 1 . List of characteristics of the orthogonally polarised modes of analysed pulsars 


Pulsar 

0144+59 

0355+54 

0450+55 

0525+21 

0809+74 

0823+26 

0834+06 

0919+06 

0950+08 

Mode-changing 

V 

V 

V 

- 

- 

- 

- 

- 

- 

Modes drift 

- 

- 

- 

- 

- 

- 

- 

V 

- 

Power law index 

of ratio of modes 

0.502 

±0.02 

0.315 

±0.004 

_ 

0.316 

±0.005 

_ 

0.137 

±0.003 

0.115 

±0.011 

0.198 

±0.010 

_ 

Are modes associated 
with components 

V 

V 

V 

_ 

V 

_ 

_ 

V 

_ 

Jump in the PA but 
no mode-switch 

_ 

_ 

V 

_ 

V 

V 

_ 

_ 

V 

Period (s) 

0.20 

0.16 

0.34 

3.75 

1.29 

0.53 

1.27 

0.43 

0.25 

Type 

T 

St 

T 

D 

Sd 

St 

D 

T 

Sd 

a (°) 

- 

51 

32 

21 

9 

84 

30 

48 

12 

yS(°) 

- 

+4.4 

+3.3 

+0.6 

+4.5 

+ 1.9 

+3.4 

+4.8 

+8.5 

w 5 o (mP) 

35 

25 

24 

48 

32 

13 

18 

23 

42 

wio (mP) 

51 

69 

82 

56 

72 

28 

25 

55 

86 

S 400 (mJy) 

6 

55 

60 

60 

80 

65 

85 

60 

400 

5 1400 (mJy) 

2 

25 

13 

9 

10 

10 

4 

4 

85 

log£(ergs/s) 

33.1 

34.7 

33.4 

31.5 

30.5 

32.7 

32.1 

33.8 

32.7 

Pulsar 

1039-19 

1133+16 

1237+25 

1737-30 

1800-21 

1929+10 

1944+17 

2016+28 

2020+28 

Mode-changing 

- 

V 

- 

- 

V 

- 

V 

V 

V 

Modes drift 

- 

- 

- 

- 

- 

- 

- 

- 

V 

Power law index 

of ratio of modes 

0.32 

±0.03 

0.234 
±0.002 [1] 

0.248 
±0.006 [1] 

2.82 

±0.05 

_ 

_ 

_ 

_ 

_ 

Are modes associated 
with components 

_ 

_ 

V 

V 

V 

V 

_ 

V 

V 

Jump in the PA but 
no mode-switch 

_ 

V 

V 

_ 

_ 

_ 

_ 

_ 

_ 

Period (s) 

1.39 

1.19 

1.38 

0.61 

0.13 

0.23 

0.44 

0.56 

0.34 

Type 

M 

D 

M 

St 

M 

T 

Tc 

Sd 

T 

a (°) 

31 

46 

53 

- 

- 

88 

19 

39 

72 

13 n 

+ 1.7 

+4.1 

0 

- 

- 

- 

+6.1 

+7.2 

+3.6 

w 50 (mP) 

43 

23 

36 

9 

97 

27 

40 

25 

37 

W10 (mP) 

50 

31 

41 

21 

321 

66 

102 

37 

49 

S 400 (mJy) 

13 

300 

110 

16.7 

20 

250 

35 

320 

110 

S 1400 (mJy) 

4 

30 

10 

6 

15 

40 

10 

30 

38 

log£(ergs/s) 

31.1 

31.9 

31.2 

34.9 

36.3 

33.6 

31.1 

31.5 

33.3 


[1] Left component only 


of the remaining profile, as in pulsar B0144+59. In the interpre¬ 
tation of the two orthogonal modes co rrespondi ng to the ordi¬ 
nary a nd extraordinary modes fearnard & Aronsll986lfPetroval 
120011) . mode-changing is likely to be a result of the spatial sep¬ 
aration of these modes due to refraction. Figures|4]through|2l1 
show that mode-changing can occur at any frequency. A com¬ 
parison between our intensity profiles of the polarisarion modes 
of pulsars B1133 +16 and B1237+25 with the results from a 
single pulse study feacker & Rankinlll98ol) suggests that OPM 
can become visible in single pulse studies even when there is 
one dominant mode thoughout the average profile. Figures 0] 
through|^ show that in many pulsars there is a great change 
in the average profile of the two polarisation modes over fre¬ 
quency, suggesting in contrast to common belief, that refraction 
can take place at all heights of emission. 


Mode-drifting 

In two cases, the longitudinal position of one of the modes rel¬ 
ative to the other can be seen to change as a function of fre¬ 
quency (see Table Q. For pulsar B2020+28, the longitudinal 
separation between the modes as a function of frequency is 
shown (see Fig. |3j(. Assuming radius to frequency mapping, 
this drift of one of the modes can be interpreted as a change 
of refraction with emission altitude. For pulsar B0919+06 the 
mode drift is less pronounced and whereas the drifting mode 
of pulsar B2020+28 moves away from the center with increas¬ 
ing frequency, the drifting mode of pulsar B0919+06 moves to¬ 
wards the center with increasing frequency. The relative drift of 
the two modes might be correlated to other observable features, 
which have to be common features of both pulsars. Table Q 
shows the common features of these two pulsars. They both 
have a period close to 0.4 s, are of type T and have a similar 
spin-down luminosity (logL) of 33.8 and 33.3 for B0919+06 
and B2020+18, respectively. 
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Spectral behavior of the ratio of the modes 

Because most data does not contain any information on the 
observed pulsar flux, the spectral index of the modes could 
only be determined relatively. Nonetheless, Figs.l22lthrouuhl24l 
show that in many cases there is a clear trend in the ratio of 
the average intensity of the two modes. The suggested power 
laws are listed in Tabled These trends imply that the spectral 
indices of the two polarisation modes often differ from each 
other. Possible explanations are that the conversion of ordinary 
modes into extraordinary modes is frequency dependent (or in 
the case of radius to frequency mapping, height dependent), 
the two modes are damped differently, or one of the modes is 
moving out of the field of view. 

Association of modes with components 

Assuming that the different modes follow different paths in 
the magnetosphere, each mode can be responsible for several 
components in the total intensity profile. In many pulsars a 
mode can, at some frequencies, indeed be associated with one 
or more components of the average intensity (see Table |TJ. 
Figuresnithrough F I I show that this can occur at any frequency. 
This association of polarisation modes with components sug¬ 
gests that the method used to separate the two polarisation 
modes is successful. 

Jump in the PA 

In several pulsars, such as B0950+08 (see Table |TJ, the PA 
jumps (sometimes gradual) 60° through 90°, while at the same 
longitudes the ellipticity angle does not change sign and the in¬ 
tensity plots of the two modes do not allow for a smooth change 
in dominance, suggesting that there is no mode-change. Yet this 
behaviour does seem to be associated with mode-changing, be¬ 
cause the ellipticity angle usually shows a similar (but smaller) 
jump and the modes are of equal strength during such a jump. A 
partial (less than 90°) or gradual jump in the PA cannot be ex¬ 
plained with two completely orthogonal modes, since the PA is 
completely determined by the strongest mode and can therefore 
only jump by 90°. In the view of the ordinary and extraordinary 
modes being 100% polarised and responsible for all the polar¬ 
ising characteristics of pulsar emission, the partial or gradual 
jump in the PA can be explained when the two polarisation 
modes are not always completely orthogonal fed war dsl2004 . 
An assumption that the modes are not 100% polarised would 
not explain these jumps. 

Correlation with type 

It is interesting that only one of the observed features, mode 
drifting, is correlated with the pulsar’s type. And even this one 
case involves only two pulsars. Also, the many pulsars we find 
containing one component, dominated by one mode at low fre¬ 
quency, and increasing in intensity with increasing frequency, 
are of all types. This suggests that the mechanism behind the 
frequency behavior of the polarisation of pulsars is similar in 
both the core as in the (multiple) conal components. 


5. Conclusions 

We have determined the average pulse profiles of two polar¬ 
isation modes for 18 pulsars at multiple frequencies, assum¬ 
ing that the modes are completely orthogonal and 100% po¬ 
larised. We find many cases where each polarisation mode can 
be wholly associated with one component in the average in¬ 
tensity and cases where a 90° phase jump in the PA is asso¬ 
ciated with the two polarisation modes changing dominance. 
This suggests that the two polarisation modes are successfully 
determined in these cases. In other cases where there is a (grad¬ 
ual) jump in the PA of less than 90° it would appear that one 
of the assumptions enabling the determination of the modes is 
not always valid. The spectra of the ratio of the integrated in¬ 
tensity of the polarisation modes show that in many cases there 
is a trend for one of the modes to become stronger with in¬ 
creasing frequency. This can be both the mode which we have 
defined as strongest, as well as the weakest mode. We also find 
that the average profiles of the modes often differ from each 
other at different frequencies. In particular, we find that when 
a component is dominated by one mode at low frequency it 
tends to increase in intensity with increasing frequency with 
respect to the rest of the profile. Furthermore, out of 18 pul¬ 
sars there are 2 pulsars that clearly show one polarisation mode 
changing its longitudinal position with frequency. From our re¬ 
sults we cannot determine whether the longitudinal transitions 
of one polarisation mode into the other at a fixed frequency are 
the result of t he X- and O-modes be ams being separated due 
to refraction iBarnard & Aronsll 198F)I) . or due to switching be¬ 
tween significant and insignificant conversion of O-mode into 
X-mode ( iPetroval200ll ). Both explanations can also account for 
some of the complex changes that occur in the average profiles 
of the polarisation modes over frequency, since plasma waves 
with different frequency will traverse different paths through 
the pulsar magnetosphere. In the model of lPetrova ( 2001 ) com¬ 
plex changes in the average profile of the polarisation modes 
over frequency implies that there are differences between the 
plasma number densities in polarisation limiting regions cor¬ 
responding to emission at different frequ encies. Following the 
technique as described in lPetroval ( 2003 ) and assuming radius 
to frequency mapping, calculation of the plasma density distri¬ 
butions from polarisation profiles at different frequencies will 
then yield changes in the plasma density both as a function of 
longitude as well as altitude. 
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Fig. 22. The ratio of the average intensity of two orthogonally polarised modes as a function of frequency, plotted on a log-log 
scale, for all pulsars. The error is calculated as 07 = Viv • cr, where N is the number of bins of the profile and cr is the noise level 
of the average intensity profile. The line through the points is the best straight line fit. Note that the scales on the vertical axis are 
chosen to be the same as to allow for comparison of the slopes. 
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Fig. 23. See Fig.[22]for explanation. Note that for PSR B1737-30 the vertical axis is shifted with respect to the other plots. This 
does not affect the slope. 
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does not affect the slope. 




















